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TECHNICAL MEMORANDUM X-53434 


ANALYSIS OF AN INTERPLANETARY TRAJECTORY TARGETING TECHNIQUE 
WITH APPLICATION TO A 1975 VENUS FLYBY MISSION 

SUMMARY 


A trajectory technique is discussed that, when systematically applied, 
enables the trajectory analyst to obtain a continuous, free-flight inte¬ 
grated trajectory from planet A to a desired target planet C via some 
intermediate target planet B. Basically, the scheme defines targeting 
parameters at the intermediate planet B in terms of the desired values 
at planet C. This allows an actual search between planets A and B, 
while, in reality, a targeting at planet C is taking place. An appli¬ 
cation of this targeting technique to a Venus flyby trajectory, i.e., 

A = Earth, B = Venus, C = Earth, is described in this report. 


I. INTRODUCTION 

The interplanetary trajectories of interest here are the flyby and 
swingby type trajectories. These trajectories, in general, depart from 
planet A, pass close by and are perturbed by an intermediate planet B, 
and arrive at a target planet C. The JPL Space Trajectories Program* 
uses a linear search routine, but a linear search between planets A and 
C is not practical because of the perturbation of the trajectory by 
planet B. Hand perturbations of the independent parameters at planet A 
are ineffective in adjusting target parameters at C. The targeting 
technique presented here will allow target parameters at planet B to 
be defined such that a targeting at C will occur. This allows the use 
of the present JPL linear search routine to be used for targeting at 
planet B. 

The targeting technique is applied to the 1975 Venus conjunction 
flyby mission where A = Earth, B = Venus, and C = Earth. 


* This JPL computer program is a high accuracy program employing either 
Encke's or Cowell's method of integration as desired. It uses precision 
ephemerldes as described by JPL Technical Release No. 34-239 and an 
oblate earth potential function containing the second, third, and fourth 
spherical harmonics. Gravitational effects of Sun, Venus, Earth, Moon, 
Mars, and Jupiter are considered. This program is described in detail 
in JPL Technical Report No. 32-223. 



II. TARGETING TECHNIQUE 


Flyby or swingby Interplanetary trajectories depart from planet A*, 
pass close by an intermediate planet B*, and arrive at a target planet C*. 

If A = C, the trajectory is called a B flyby trajectory, and if A / C, 

the trajectory is called a B swingby trajectory. Figure 1 shows a typical 
Venus flyby trajectory where A “ Earth, B = Venus, and C = Earth. 

Knowing the relation between the planets of interest (dates, helio¬ 
centric longitude, etc.) and the approximate trip times from A to B and 
B to C, we can easily find the flyby or swingby trajectories by two-body 
approximating conic programs [6, 7]. Using the conic program, we deter¬ 
mine the flyby or swingby trajectories by analyzing groups of A to B 
trajectories and B to C trajectories. The conditions necessary at B 
to insure a flyby or swingby trajectory are as follows: (1) the time 
of arrival at B from the A to B trajectory must equal the depart time 
for the B to C trajectory and (2) the arrival hyperbolic excess velocity 

magnitude from the A to B trajectory at B must equal the depart hyper¬ 

bolic excess velocity magnitude for the B to C trajectory. With these 
conditions met, the angle between the arrival hyperbolic excess velocity 
vector and the depart hyperbolic excess velocity vector, which is the 
bend angle, determines the radius of closest approach (RCA) at B. With 
conditions (1) and (2) met at B, the trajectory from A to B is said to 
be merged with the trajectory from B to C, thus yielding a flyby or 
swingby trajectory. Reference 6 describes a program which has the above 
conditions automated, and the selection of RCA is left for the analyst. 
Venus swingby trajectory opportunities with respect to Earth-Mars and 
Mars-Earth trajectories can be found in References 2 and 9 and flyby 
trajectories to Mars and Venus can be found by using References 8 and 10. 

TWo-body approxiraat;Lng conic trajectories are good for preliminary 
analysis. However, for detail mission analysis where communication 
distances, communication angles, tracking requirements, and point by 
point time histories are needed, a precision integrated trajectory 
becomes necessary. Here, in particular, the nominal or reference free 
flight trajectory is of importance. 

The JPL Space Trajectory Program [3, 5] is used to determine the 
integrated flyby or swingby trajectory. As presently used, the JPL pro¬ 
gram employs a targeting method which provides the capability of finding 
A to B type trajectories. To clarify terminology, the JPL targeting 
method will be discussed briefly. 


* A, B, and C are used to refer to planet A, planet B, and planet C, 
respectively. 
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The JPL trajectory targeting method is the N x N Newton-Raphson 
scheme as defined in Reference 5. In this scheme, a nominal trajectory 
to B, usually using conditions from approximating conic trajectories, 
is made to establish a column matrix of miss components, P. The 
partials are obtained by N perturbed trajectories and differencing. 
Having obtained a matrix M of partials, the search routine solves the 
equation 


M AX = P (1) 

to obtain incremental conditions, as represented by the column matrix AX, 
at A to null the miss components at B. The partials found by differencing 
are used for three iterations and then are recomputed. This procedure is 
repeated until the miss components are driven to zero or within some 
convergence tolerance. 

The miss components are the differences between the desired dependent 
variables at B and the values obtained in the nominal trajectory. The 
dependent variables are usually the two components of the impact parameter 
B and the time of flight, Tp, from A to B or hyperbolic excess velocity, 

Vh» on arrival at B. The two components of the imgact parameter B are 
known as the JPL targeting system [5] or the B*T, B*R system. 

The B*T, B*R target system as used in the JPL Space Trajectory Program 
[3] will be described. In Figures 2 and 3, S is a unit vector parallel to 
the incoming asymptote referenced to the center of target or planet of 
interest, T is a unit vector parallel to or lying in the plane of interest 
(ecliptic plane, earth's equatorial plane, etc.) and perpendicular to S, 
and R completes the orthogonal system R, S, and f. The vector B lies in 
the R - T plane and perpendicular to the approach asjnmptote (see Figures 
2 and 4). B has the magnitude of the semi-major axis, b, of the approach 
hyperbola. In Figure 5, looking at the S - ^ plane, we see that 5*6 and 
and B»f are the projections of B on S and f, respectively. 

Injection conditions at A are usually referred to as the independent 
variables for the search procedure. These independent variables as used 
here are (1) T^, launch time of day which, in conjunction with a specified 
earth launch azimuth, establishes a specific departure orbit, (2) Tqq, 
the departure point or coasting time in the specified earth orbit, and 
(3) Tgu, the thrusting time of the earth escape maneuver. 

Now using the independent variables T^* Tqq, Tg^ and the dependent 
variables B*T, B*R, Vfl, the elements in equation (1) become 
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^B-f 

SB-f 

^B*T 


^'^CO 

^"^BU 

dB-R 

Sb-r 

3b*r 


^”^00 

^'^BU 

avH 

aVH 



^'^CO 

^'^BU 


M = 


ATt 


AT, 


CO 


AT- 


BU 


at planet A 


(3) 


P = 


^♦T 


AV, 


H 


at planet B 


(4) 


The procedure for obtaining these values is as described above. 

Since the flyby or swingby trajectory is hyperbolic about B, the 
set of dependent variables above are linear for reasonable perturbations 
of the independent variables. This linearity arises from the fact that, 
with small perturbations of the independent variables at A, the hyperbolic 
excess velocity vector is translated and not rotated at planet B [5J. 
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After finding the desired conditions at planet B, as specified from 
the approximating conic conditions, the trajectory is allowed to continue 
by B, and its relationship to planet C is observed. From Reference 10 
and from experience in obtaining flyby trajectories, the first try usually 
has a large RCA at C (approximately 6 x 10® km for Venus flyby and Mars 
flyby trajectories where C = Earth). Attempts to move the conditions at 
B to reduce the RCA at C by using various curve-fitting schemes and manual 
perturbations of independent variables at A yield trajectories which 
return to C with first order of magnitude improvements (3 x 10® km for 
Venus flyby and Mars flyby trajectories where C = Earth). 

Reference 4 has shown that the effects at C from errors occurring 
in hyperbolic excess velocity (Vg), right ascension (a) and declination 
( 6 ) of the incoming asymptote, and time of flight from A to B (T^-g) at 
B_are small in comparison to errors in the impact parameter components 
(B*T and B*R) at B. Targeting at C from flyby or swingby conditions at 
B will be described in terms of the impact parameter components at B and 
C. The effects of small perturbations in the impact parameter components 
at B on the impact parameter component at C must be found and used in a 
manner to obtain desired changes at B to null the errors in the impact 
parameter components at C, thus targeting the trajectory at C. From 
Reference 4, the variation in impact parameter components at B and their 
effects on target parameters at C are given by the R matrix: 


as-Tc 

as-fc 

as-Tg 

as-Rg 

ai-Rc 

aS-Rc 

as.Tg 

aS-Rg 


(5) 


By knowing the R matrix, changes in impact parameter components at 
B produce associated changes in impact parameter components at C. 


^•Tc 


R 


Z^-Rg 


( 6 ) 


** Subscripts B and C indicate planet B and C, respectively. 
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However, the desired conditions at C are known and the conditions at B 
are desired in order to target the flyby or swingby trajectory at B. 




• ^ 

<!^*Tb 


^•Tc 


II 

I 

H 


/N 


— /\ 

Z!B*Rb 


AB-Rc 


Using equations (7) and (1), and the R matrix, we can find the 
desired changes in impact parameter components at C, and finally the 
needed changes in impact components at B. The process for finding the 
"best" trajectory is an iterative process. "Best" is used to mean 
(1) convergence to the desired conditions at C or (2) that the desired 
changes at B fall within a convergence tolerance. 

Figure 6 shows a simplified logic flow diagram of the above procedure 
for obtaining targeting at C using the JPL Space Trajectory Program. 

In the next section, the targeting technique is used to obtain a 
flyby trajectory where A = Earth, B = Venus and C = Earth. Only the 
targeting technique from B to C will be discussed in detail, since 
adequate documentation of the JPL targeting method is found in References 
3, 5, and 10. 


III. VENUS FLYBY TRAJECTORY 

In phase I of the Manned Venus/Mars Flyby Study [10], minimum 
weight in earth-orbit missions for various Venus conjunctions are 
presented. The minimum weight in Earth orbit mission for the 1975 
Venus conjunction was selected from Reference 10 to illustrate the 
targeting technique where A = Earth, B = Venus, and C = Earth and to 
provide a reference trajectory for phase II of the Manned Venus/Mars 
Flyby Study for the 1975 Venus conjunction. 

Using the two-body conic program to establish a preliminary trajec¬ 
tory, a trajectory having the characteristics given in Table 1 was chosen. 
This is a free flight trajectory which departs Earth (A) on June 7, 1975, 
passes Venus (B) on October 2, 1975, and returns to Earth (C) on June 7, 
1976, for a total trip time of 366.4 days. A typical 1975 Venus trajec¬ 
tory is shown in Figure 1. 
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The conic arrival and depart conditions at Venus were used to estab¬ 
lish the B’R^and for the initial dependent search parameters. 

The JPL trajectory targeting method was used with the independent search 


parameters, T^, 


Tco> 3.nd T 3 U. 


A converged trajectory from Earth was 


obtained. This trajectory was allowed to continue past Venus and to 
return toward the Earth. The RCA at Earth return was 5,696,092 km. 
Conditions on this first pass trajectory are given in Table 2, Column 1. 


The perturbed trajectories were found using the same procedure as 
discussed above end are tabulated in Table 2, Columns 2 and 3. With 
these conditions, the R matrix was computed and is tabulated in Table 2B 
as First R Matrix. Using this R matrix and equation (7) of Section II, 
the and z5*&^were calculated and are shown in Table 2, Column 4. 

(Note: Desired B*T® = 13430 km and B«R* = 50 km.) With a new B»T^ and 



a trajectory was found as shown in Table 2, 


Column 5 where the RCA 


at Earth is 687,691 km. This procedure continued for two more iterations 
as seen in Table 2. A new R matrix was then calculated (Table 2B, Second 
R Matrix) and used for one iteration (Table 2A, Column 5). This trajec¬ 
tory, Table 2A, Column 5, yielded a and Table 2A, Column 6, 


less than the convergence tolerance of 1.5 km used in the JPL Space 
Trajectories Program [3]. The "best" trajectory was then printed in 
detail. 


The conic trajectory and integrated trajectory are compared in 
Table 1. The hyperbolic excess velocities leaving Earth, arriving and 
departing Venus, and arriving Earth are in close agreement. This agree¬ 
ment confirms that conic approximation assumptions used in calculating 
weight in Earth orbit of Reference 10 are valid. 

For purposes of communication between the spacecraft and Earth, for 
the 1975 Venus flyby trajectory. Figures 7, 8, and 9 show the pertinent 
communication angle and distances from Earth and Venus. A maximum 
distance from Earth of 117,360,440 km occurs 190 days after injection. 

Table 3 is a tabulation of the thrusting trajectory for the Earth 
escape maneuver. The initial weight is based on Reference 10 where two 
J-2 engines are used in the Earth escape stage for the minimum weight 
Venus mission during the 1975 conjunction. Thrusting is along the 
vehicle's longitudinal axis and in the plane of motion. This trajectory 
starts from a 485 km orbit on June 7, 1975 at 6 hours, 49 minutes, and 
27.591 seconds past midnight, GMT, The 485 km Earth orbit was chosen 


* T is referenced to Earth equatorial plane. Subscripts ^ and ® refer 
to Venus and Earth, respectively. 
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-■0 * 


for its rendezvous and lifetime characteristics, both of which are 
important for manned missions involving earth orbital operation. 

Figure 10 illustrates the coordinate systems used in the trajectory 
tabulations. 

Table 4 is a listing of various trajectory parameters for the Venus 
flyby trajectory. The trajectory is tabulated in five phases: (1) Earth 
Depart, (2) Heliocentric Earth-Venus, (3) Aphrodiocentric, (4) Heliocentric 
Venus-Earth, and (5) Earth Return. After a thrusting period of 393.6 
seconds (Table 3), the spacecraft reaches mission injection on June 7, 

1965, at 6 hours, 56 minutes, and 1.192 seconds, GMT, with a hyperbolic 
excess velocity of .1091 EMOS. The time from injection to pericenter 
passage at Venus is 116 days, 0 hours, 34 minutes, and 29.06 seconds. 

The trajectory passes Venus with an RCA of 6354.2 km on October 1, 
1975, at 7 hours, 30 minutes, and 30.252 seconds, GMT, and returns to a 
geocentric radius of 9717.2 km with a hyperbolic excess velocity of 
.2544 EMOS on June 7, 1976 at 4 hours, 58 minutes, and 27.824 seconds, 

GMT. A total mission trip time of 365 days, 22 hours, 2 minutes, and 
26.6 seconds is required. 


CONCLUSIONS 

The targeting technique presented here provides a means whereby 
accurate Integrated reference trajectories for flyby or swingby missions 
can be obtained. For the Venus swingby illustration , the RCA on arrival 
at Earth was improved by an order of magnitude over that obtained using 
previous techniques. 

IBM 7094 computer time for obtaining the Venus trajectory using this 
targeting technique was approximately 4.1 hours. The R matrix and incre¬ 
ment values were calculated by hand. This, however, reduced computation 
time by a factor of four over that required for previous techniques. 
Therefore, the targeting method described in this report allows one to 
obtain a more desirable trajectory with a considerable reduction in the 
amount of computer time. 
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TYPICAL VENUS FLYBY MISSION PROFILE FOR 1975 CONJUNCTION 
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FIG. 4. B IMPACT PARAMETER 
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AND B-R IN THE R-T PLANE 
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Flight Timt (doys) 

FIG. 7. VENUS FLYBY MISSION 1975 CONJUNCTION 
VEHICLE-EARTH ANGLE AND SUN-EARTH-VEHICLE ANGLE 

VERSUS FLIGHT TIME 








Venus Possoge 
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FIG. 8. VENUS FLYBY MISSION 1975 CONJUNCTION 
SUN-VEHICLE-VENUS ANGLE AND SUN-VENUS-VEHICLE ANGLE 

VERSUS FLIGHT TIME 






Distance (AU) 
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FIG. 9. VENUS FLYBY MISSION 1975 CONJUNCTION 
EARTH-VEHICLE DISTANCE AND VENUS-VEHICLE DISTANCE 

VERSUS FLIGHT TIME 







Venus Orbit Plone 



18 


FIG. 10. COORDINATE SYSTEMS 











TABLE 1 


CCMPARISON BETWEEN INTEGRATED AND CONIC CALCULATIONS 


V 

CONIC 

INTEGRATED 

Launch Date (J. D.) 

2442570.0 

2442570.5 
(Jun 7, 1975) 

Earth-Venus Trip Time (days) 

117.0 

116.0 

Hyperbolic Excess Velocity 

Leaving Earth (EMOS)* 

.1086 

.1091 

Right Ascension of Incoming 

Asymptote at Venus (deg) 

180.140 

180.510 

Declination of Incoming 

As 3 nnptote at Venus (deg) 

3-771 

4.287 

Hyperbolic Excess Velocity 

Approaching and Leaving Venus (EMOS) 

0.1554 

0.1554 

Right Ascension of Outgoing 

As 3 niiptote at Venus (deg) 

91.448 

89.042 

Declination of Outgoing 

Asymptote at Venus (deg) 

22.193 

24.056 

Venus-Earth Trip Time (days) 

249.4 

249.9 

Hyperbolic Excess Velocity 

Arriving Earth (EMOS) 

0.2546 

0.2544 


NOTE: All angles listed above are referenced to aphrodiocentric 
coordinate system with Xv-Yv plane parallel to the earth's 
equatorial plane 

*EMOS = Earth Mean Orbital Speed (29.7849 km/sec) 
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Tabulated Data Using First R Matrix 
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TABLE 2A 

Tabulated Data Using Second R Matrix 
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TABLE 2B 


Listing of R Matrix Elements 

First R Matrix 



B . 

B . R, 

B • T« 

+9218.57 

-4871.73 

B • 

- 245.93 

+ 992.90 


Second R Matrix 


a 

B ♦ 

1 

B • Rg 

B • 

+9652 

.26 

-3213.41 

B • 

+ 189 

.39 

862.79 
























TABLE 3A 


THRUSTING TRAJECTORY FOR EARTH ESCAPE 


Definition of tabulated trajectory values contained in the thrusting 
trajectory for Earth escape. 


T time from June 7, 1975 (6 hours, 49 minutes, 27.591 sec) 


Line #1 


Geocentric Equatorial 



Vehicle Cartesian position components. (Coor¬ 
dinate system origin at center of earth, X-axis 
in direction of earth's of date vernal equinox, 
X-Y plane is in the earth's equatorial plane) 
(km) 



Vehicle Cartesian velocity components measured 
in same coordinate system as above (km/sec) 


Line #2 Geocentric Equatorial 

RE Radius from center of earth to vehicle (km) 

LATE Geocentric latitude (deg) 

LONGE Geocentric longitude (deg) 

VE Earth-fixed speed of vehicle (km/sec) 

PTE Earth-fixed path angle from the horizontal (deg) 

AZE Earth-fixed azimuth angle from north (deg) 
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END OF INJECTION THRUSTING _ BURN OUT WEIGHT 264,652 lbs. (120,044 kg) 

INJECTED PAYLOAD 172, 583 lbs (78, 282 kg.) 









TABLE 4A 


1975 VENUS-JLYBY TRAJECTORY 
Geocentric 


Definition of tabulated trajectory values contained In the Earth 
Depart trajectory and Earth Return trajectory. 

T = time from injection from earth (days) 

Line #1 Geocentric Equatorial 


Xe 

Ye 

Ze 


} 


Vehicle Cartesian position components. (Coordinate 
system origin at center of earth, X-axis in direction 
of earth's of-date vernal equinox, X-Y plane of 
earth's equator) (km) 


Xe 

Ye 

Ze 


) 


Vehicle Cartesian velocity components measured 
in same coordinate system as above (km/sec) 


Line #2 Geocentric Equatorial 

Re Radius from center of earth to spacecraft (km) 

LATe Geocentric Latitude (deg) 

LONGe Geocentric Longitude (deg) 

Ve Earth-fixed speed of spacecraft (km/sec) 

PTe Earth-fixed path angle from the horizontal (deg) 

AZe Earth-fixed azimuth angle from north (deg) 
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Line #3 Heliocentric Ecliptic* 


Rs Radius from center of sun of spacecraft (km) 

LATs Celestial latitude of spacecraft (deg) 

LONGS Celestial longitude of spacecraft (deg) 

Vs Inertial speed of spacecraft (km/sec) 

PTs Path angle from the normal to the heliocentric 

radius vector (deg) 

AZs Azimuth angle from celestial north (deg) 

Heliocentric 


Definition of tabulated trajectory values contained in the helio¬ 
centric Earth-Venus transfer trajectory and heliocentric Venus-Earth 
transfer trajectory. 


T = time from injection from earth (days) 
Line #1 Heliocentric Ecliptic 


Xs 

Ys 

Zs 


Vehicle Cartesian position components. (Coordinate 
system origin at center of sun, X-axis in direction 
of earth's of-date vernal equinox, X-Y plane in the 
ecliptic) (km) 


Xs 

Ys 

Zs 


} 


Cartesian velocity components measured in same 
coordinate system as above (km/sec) 


Line #2 (Same as Geocentric Line #2) 


Line #3 (Same as Geocentric Line #3) 


* (Origin of coordinate system at center of sun, X-axis 
earth's of-date vernal equinox, X-Y plane in ecliptic 


in direction of 
plane.) 
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Aphrodlocentrlc 


Definition of tabulated trajectory values contained in the Venus 
passage trajectory. 

T = time from injection from earth (days) 

Line #1 Aphrodiocentric Earth Equatorial 


Xv 

Yv 

Zv 


) 


Vehicle Cartesian position components. (Coordinate 
system origin at center of Venus, X-axis in the 
direction of earth's of-date vernal equinox, X-Y 
plane of earth's equator) (km) 


Xv 

Yv 

Xv 


Vehicle Cartesian velocity components measured 
in same coordinate system as above (km/sec) 


Line #2 (Same as Geocentric Line #2) 

Line #3 (Same as Geocentric Line #3) 


Line #4 Aphrodiocentric Earth Equatorial* 

Rv Radius from center of Venus to spacecraft (km) 

DECv Declination (deg) 

RAv Right ascension (deg) 

Vv Speed relative to Venus (km/sec) 

PTv Path angle from the horizontal (deg) 

AZv Azimuth angle as referenced to geocentric equatorial 

coordinate system (deg) 


* Reference plane is the earth's equator plane. 
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1975 VENUS FLYBY TRAJECTORY 
EARTH DEPART TRAJECTORY 
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* PERI-CENTER PASSAGE AT VENUS ON OCTOBER 1, 1975 07 HRS 30 MIN 30.252 SEC. (GMT.) 
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